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ABSTRACT 

We study the possibility of spontaneous CP violation (SCPV) at 
the tree level in models with an extended Higgs sector. We show 
that the minimum equations for the complex phases of the vacuum 
expectation values (VEVs) have always a geometrical interpreta- 
tion in terms of triangles. To illustrate our method we analyze 
the minimal supersymmetric (SUSY) model with R-parity violat- 
ing couplings and sneutrino VEVs, where there is no SCPV. Then 
we study SUSY models with extra Higgs doublets and/or gauge sin- 
glets, and find that the simplest scenario with SCPV must include 
at least two singlet fields. 



1 Introduction 



Although the observed CP violation can easily be accommodated in the standard model, 
allowing for arbitrary CP violating terms leads to phenomeno logical difficulties in models with 
an enlarged Higgs sector. For example, in general two-Higgs doublet models with acceptable 
flavor changing interactions the prediction for e# would still be too large Q . Another example 
is the minimal SUSY extension of the standard model (MSSM), where arbitrary complex 
phases in gaugino masses and scalar trilinears would produce too large electric dipole moments 

In order to bring the predicted CP violation in such extensions to acceptable levels one 
can consider imposing additional symmetries, like a flavor symmetry or a discrete symmetry 
in the Higgs sector. One obvious possibility is to impose CP invariance itself. If CP violation 
is spontaneous ||, in the sense that it appears via VEVs of scalar fields, then some phases 
could be naturally suppressed by the ratio of two mass scales present in the model, like the 
top and bottom quark masses 0, |5[] or the electroweak and the SUSY scales. This hierarchy 
could be used, for example, to accommodate small complex phases in gaugino masses and 
scalar trilinears together with the large CKM phase required to explain the kaon system. 

Another interesting motivation for SCPV arises from the strong CP problem. As far 
as CP is a good symmetry the QCD phase 9 is zero. The spontaneous breaking of CP could 
respect this initial value of 9 while creating the observed weak CP violation || [7]. Such 
scenarios find a natural framework in left-right models, where it was recently shown that a 
minimal model of this type must be supersymmetric and with a low scale of SU(2)r symmetry 
breaking |§. 

In consequence, SCPV appears as a well motivated possibility in SUSY models. A 
priori, these models contain enough ingredients and arbitrarity to introduce SCPV: a neutral 
scalar sector with at least two Higgs fields plus three sneutrinos, and a large number of 
arbitrary SUSY-breaking terms. However, it is well known that the MSSM has only real 
minima. In fact, complex VEVs are possible once radiative corrections are included ||, 
but then the model contains a too light Higgs boson which is experimentally excluded [TO 



Therefore, it is necessary to introduce other fields and couplings. Two cases have been 
considered before: an extra gauge singlet field (2D IS model) and an extra pair of Higgs 
doublets (4D model). The presence of a singlet is appealing because it defines a scenario 
where the SUSY mass term of the Higgs is substituted by a trilinear term in the superpotential 
combined with a singlet VEV, avoiding the n problem and relaxing the upper bound on the 



mass of the lightest neutral Higgs [O]. The singlet model (known as the next-to-MSSM) has 



been analyzed in this context by Romao ||12|| , with the result that no CP violating minima 
exist at tree levelQ. The inclusion of an extra pair of Higgs doublets is also an obvious 
generalization of the MSSM. But the result there is negative as well, the minimum in the 4D 



model is always real |H |. In Section 4 we review the arguments used to prove these results. 

Contrary to non-SUSY models, where the Higgs quartic interactions are arbitrary, 
SCPV in SUSY scenarios seems to require an increasing number of species and fields. This 
fact quickly will make the analysis of minimum equations cumbersome. In Section 2 we 
describe a method to find complex minima that interprets the equations for the phases in 
terms of a geometrical object combining triangles. The method is general in the sense that it 
can be applied to any potential (SUSY or not). We illustrate it analyzing a SUSY model with 
minimal matter content but R-parity breaking couplings and sneutrino VEVs. Our objective 
in this paper is to find the minimal SUSY scenario giving SCPV in the Higgs sector at the 
tree level. In Section 3, we extend the Higgs sector of the MSSM with singlets and/or extra 
doublets (the addition of SU{2)l triplets is phenomenologically disfavored), and we show 
that the minimal scenario consists of at least two singlets, regardless of the number of Higgs 
doublets. In Section 4 we give our final remarks and conclusions. 



2 The geometrical method 

In this section we describe a geometrical method of analyzing the minima equations for the 
phases that was recently used in the context of a SUSY model with four Higgs doublets. 
The method is actually a generalization of the procedure used in the simplest case, where 



the object representing the equations is just a triangle O, 15]. Let us start studying the 



non-SUSY model proposed in [16] to review this simplest case, and then we will describe in 



some detail a more general scenario. 



Consider the extension of the standard model with three Higgs doublets [IJJ where 
one doublet couples to the up quarks, another one to the down quarks, and the third one 
does not couple to quarks at all. The scalar potential is given by 

V 3D = m 2 i H}H i + KiHjH,) 2 + X'^HjH^^H,) 
+ [\ 12 {H\H 2 ){H\H 2 ) + \ 13 {H\H 3 ){H\H 3 ) 
+\ 23 {H\H 3 ){HlH 3 ) + h.c] , (1) 



: We do not consider here singlet extensions where dimcnsionful couplings are allowed. In such case it was 
shown that spontaneous CP violation was possible |f3| . 



with all the parameters real. The neutral components of the doublets will have VEVs 

(tf^-L^e* (z = l,2,3), (2) 

where an hypercharge transformation is used to set S\ = (the minimum will be degenerate 
due to the U(l)y symmetry of the potential). The value of vi and 8i will be given by the 
solution of the minimum equations. For the complex phases these are (||^ = 0) 

\i 2 vfv% sin 8 2 + \23V2 V 3 sin(<5 2 + 8 3 ) = 

X 23 vlvl sin(5 2 + 8 3 ) + \13vfvl sin 83 = (3) 



To solve these equations we draw a triangle |T5[ with sides a x , a 2 , a 3 and opposite angles 
7r — 8 2 , 7r — 83, 8 2 + 83 — n, respectively (this is the lower right triangle in Fig. 1). If we define 

0-3 = A 2 3f2^3 ( 4 ) 

then the sine law applied to the triangle implies Eq. (3), i.e., the triangle is the solution to 
Eq. (3). Using the triangle it is now possible to give the values of the cosines of the phases 
(which appear in the minimum equations for the moduli) in terms of the sides: 



A 12 t>it>2 ^32vlvl Af 3 f?f|' 

r (A 32 ^ 2 2 )(Ai 2 ^ 2 2 ) 1 2 I a 1 , 

3 ^ 2 K A 13 ^r l A 32 ^ 2 2J [ %2vlvl 

(Xl2vlv^)(X 13 vlvl) 1 2 1 2 1 , 



cos(<5 2 + 83) 



\32V3Vl \\2v\v\ A? 3 «?V|' 



(5) 



Substituting these expressions in the three minimum equations for the moduli we will obtain 
the equations in terms of the three moduli i>, only (with no phases), and these equations 
can be solved numerically. The value of the moduli will fix the sides of the triangle and, in 
consequence, the value of the complex phases. For the particular three doublet model under 
study, it is easy to find solutions for adequate values of the couplings, and then there is 
SCPV. Below we analyze a similar example where the phases can also be expressed in terms 
of the moduli, but then the minimum equations for the moduli are always incompatible. 



This method of expressing the phases as functions of the moduli can be generalized to 
models with more than two complex phases. We will show now how to build a set of triangles, 



with their angles related to the phases in the VEVs, that solves the phase minimum equations. 
The method can be used if the equations involve only sines (and not cosines) of the phases, 
which is the case for Higgs potentials with all the couplings real. In models with CP odd 
fields the CP invariant potential may include purely imaginary couplings. However, since in 
SUSY models all the scalar fields are complex, the factors of i in the couplings can always be 
rotated away by field redefinitions. In non-SUSY cases there may be exceptions with CP odd 
real scalarsQ. The (real) VEVs of such fields would break CP, which would be transmitted to 
the rest of the Lagrangian through complex couplings either through additional fermions 



or through explicit couplings to the Higgs doublets. The minimal model of the latter type is 
a two Higgs doublet model with flavor conserving couplings [19j and with one CP odd real 
singlet. In this case the minimum of the potential has terms sin(52) and cos(2#2), being <5 2 
the relative phase of the two doublets. Here, however, we can also do a rotation of 82 by 7r/2 
to convert both terms into cosines, and the minimum equation will only involve sines of the 
phase. 

To describe our method in some detail we will consider an extension of the MSSM 
with R parity violating couplings (R model). The model has a minimal matter content with 
one pair of Higgs doublets. However, sneutrino VEVs are allowed. This fact will introduce, 
in an effective way, three new doublets of —1/2 hypercharge in the scalar Higgs sector. The 
relevant part of the superpotential is 

W = W M ssm + Hi H 2 U (6) 

We will restrict our analysis to the VEVs of the Higgs (Hi, H 2 ) and lepton (Lj) fields, assuming 
that we are in the charge conserving part of the parameter space. We will closely follow the 
notation used in the multi-Higgs case in fBfl , with odd (even) indices indicating —1/2 (+1/2) 
hypercharge. We redefine (Li,L 2 ,L 3 ) — > (H 3 ,H 5 ,H 7 ) and write the VEVs for the neutral 
components of these doublets as 

(Hi) = ±=v^ (1 = 1,2, 3, 5, 7), (7) 

using a global hypercharge transformation to set 5i — 0. Then the VEV of the scalar potential 
is 

t t 1 2 2 1 1,2 2 1 1 2 2 1 1 2 2 , 1 2 2 , 2 r 

V V = 2 m l V l + ^ m 2 v 2 + 2 m 3 V 3 + 2 m 5 V 5 + 2 m 7 V 7 + m 12^2 COS 2 

+ml 3 v iv 3 cos #3 + ml 5 viV5 cos £5 + m\ 7 viVj cos 5-j + rr? 32 v z v 2 cos (£3 + 82) 
+ml 2 v 5 v 2 cos(<5 5 + 5 2 ) + rrij 2 v 7 V2 cos(5 7 + 5 2 ) + ml 5 v 3 v 5 cos(<5 3 - 5 5 ) 
+ml 7 v 3 v 7 cos(5 3 - 5 7 ) + ml 7 v 5 v 7 cos(5 5 - 5 7 ) + V D 



2 We thank Goran Senjanovic for pointing out this posssibility. 



where all the mass parameters are real and Vd is the D-term part of the potential 

Vd = ^(9 2 + 9' 2 ) H + v! + v\ + v 2 - vlf (9) 



We are seeking CP- violating minima, i.e., minima where some phases are different 

av 

c)5, 



from or tc. First, we solve the minimum conditions for the phases (—It: = 0) 



m\ 2 viv 2 sin S 2 + ml 2 v 3 v 2 sin(<5 3 + 5 2 ) + m\ 2 v 5 v 2 sin(<5 5 + 5 2 ) + m 72 v 7 v 2 sin(<5 7 + 8 2 ) = 
m 2 32 v 3 v 2 sin(<5 3 + 5 2 ) + mj^vs sin 5 3 + m 2 35 v 3 v 5 sin(<5 3 - 5 5 ) + m 2 37 v 3 v 7 sin(<5 3 - 8 7 ) = 
^52^5^2 sin(<5 5 + 5 2 ) - ml 5 v 3 v 5 sin(<5 3 - 5 5 ) + m 2 lb viv 5 sin 5 5 + m 2 b7 v^v 7 sin(<5 5 - 5 7 ) = 
m 2 2 v 7 v 2 sin(5 7 + 5 2 ) — m\ 7 v 3 vj sin(5 3 — 5j) — m^ 7 v^v-j sin(5 5 — 5 7 ) + m\ 7 viv 7 sin 5 7 = 

(10) 

Now we express the solution of the minimum equations for the phases in terms of a com- 
bination of triangles. We observe that there are ten independent quantities m^ViVj in the 
equations and, in consequence, the space of solutions will be 10-dimensional. The geometrical 
solution must be given by a set of six triangles whose angles involve the different combinations 
of four phases in Eq. (0). Such a combination of triangles contains ten independent sidesQ. 
A choice for the six triangles is shown in Fig. 1. Adding the sine law applied to the triangles 
it is straightforward to obtain the set of equations in (p!U|), where we identify 

m\ 2 viv 2 = ai + hi + x% , m\ 2 v 3 v 2 = a 3 , m 2 52 v 5 v 2 = -b 3 , 

m 2 l3 viv 3 = a 2 + ci + z x , m 2 2 v 7 v 2 = -x 3 , m 2 35 v 3 v 5 = -c 3 , 

^15^5 = h - c 2 + yi , m 2 37 v 3 v 7 = -z 3 , ml 7 v 5 v 7 = -y 3 . 

ml 7 v x v 7 = x 2 - z 2 - y 2 , (11) 

Now we have to express cos5j in terms of vn^ViVj and substitute them in the five minimum 
equations for the moduli. If there we find a solution fixing Vi, then we have SCPV; if the 
equations are incompatible then any minimum that may exist will be CP conserving. Al- 
though this procedure sounds simple, in general it is not easy to obtain analytic solutions. In 
particular, to express the sides of the triangles in terms of rriijViVj involves solving a quartic 
equation. The multitriangle is useful to generate numerical solutions, as we do in Section 3 
to study other models, but it is not enough to solve analytically complicated cases. In the $ 

3 Note that the four angles and one side of each triangle also fixes the six triangles. In general, the number 
of phases in the equations plus the number of triangles in the multitriangle must be equal to the number of 
independent parameters in the equations. 



model under study here, however, we can find a way to simplify the equations that allows for 
analytical solutions. We can rotate the original fields (Hi,H 3 ,H 5 ,H 7 ) to a new basis so that 



2 2 2 2 2 2 n 

m 13 = m 15 = m 17 = m 35 = m 37 = m 57 = • ( 12 ) 

This rotation redefines the diagonal terms mf but does not introduce any complex phase in 
the Lagrangian. From Eq. fli~0| ) this set of zero masses forces 



sin 5 2 = sin(5 2 + 5 3 ) = sm(5 2 + 5 5 ) = sm(5 2 + 5 7 ) = (13) 

which, if the five VEVs Vi are non-zero, forces all the phases to be trivial. Then there is no 
CP violation in the general case with all v , 7^ 0. 

If the modulus v 7 = and all the other VEVs nonzero, then the phase S 7 is irrelevant 
(from (|TT| ) we see that X3 = = Z3 = and the corresponding triangles become infinite). 
Now the three-phase case (62,63,55) that results has three triangles left and is completely 



analogous to the 4D model discussed in [14], with no SCPV. 



The case with two VEVs zero, v 5 = v 7 = 0, appeared in f20fl , where it is claimed that 
there are CP- violating minima. Since we disagree with this result, let us consider it in a bit 
more detail. Now the minimum equations for the phases reduce to 

m\ 2 viV2 sin 6 2 + ml 2 v 3 v 2 sin(5 3 + 6 2 ) = 

m\ 2 v 3 v 2 sin(5 3 + 6 2 ) + mj 3 viv 3 sin 63 = (14) 

Again, the lower right triangle in Fig. 1 solves these equations if a x = mi 2 viv 2 , a 2 = 
^13^1^3 and a 3 = m 2 3V 2 V3. The cosines of the angles can be obtained from the expressions in 
just exchanging Xijvfvf — > m^ViVj. 

Notice that this does not mean yet that there is SCPV, we still have to solve the three 
minimum equations for the moduli (i>tfj: = 0): 

m\v\ + m\ 2 viv 2 cos 6 2 + m1 3 viV 3 cos 5 3 + v\g{y) = 
m\v 2 + m 12 v iv 2 cos 5 2 + m\ 3 v 2 V3 cos(5 2 + S3) — ^^(v) = 
mlvl + 771^3^1^3 cos 63 + ml 3 v 2 V3 cos(<5 2 + 5 3 ) + vlg(\) , = 

(15) 

where g(v) = \(g 2 + g' 2 )(v\ + v 3 — f|). We substitute the expressions for cos5j in the equations 
above and obtain 

2 2 

,2r™2 m l2 m l3 



v 2 [ml _._^pZ +g{v)] = Q 



<m 2_m^s_ g(y)] = Q 

,a 12 

j M _rn$^ +g(y)] = Q 

m 13 

(16) 

There is only one combination of VEVs, g(v), to solve the system of three equations. Then, 
unless the masses are fine tuned in such a way that two of the moduli equations are trivial, the 
three equations can not be satisfied simultaneously. This implies that the minimum equations 
for the phases must have the trivial solution (phases equal zero or tc). We conclude that there 
is no tree-level SCPV in the extension of the MSSM with R-parity violating couplings. 

In the two simple non-SUSY and SUSY cases above the phase minimum equations 
have been solved in terms of triangles. It is now possible to see how to extend the procedure 
to more complicated cases. If we add a Higgs field to any of the previous models, this field 
will introduce one more phase and will increase the number of independent parameters in 
the equations. To find the solution now we need to draw new triangles, and the number 
of independent distances in these triangles must be exactly equal to the number of new 
parameters (couplings) in the equations. This will be always possible because we can draw 
an arbitrary number of new triangles with the existing phases, and each new triangle will 
introduce one more independent distance: the overall scale of the triangle, which is not fixed 
by the phases. We can always add triangles until matching the number of new couplings in 
the minimum equations. The procedure will become more obvious in the next Section, where 
we consider SUSY models with additional Higgs fields. 



3 Search for the minimal SUSY model 

The method described in the previous section has been already applied to minimal SUSY 
scenarios, namely, to the 2D1S model and to the 4D model. In both cases the results are 
negative. In the 2D1S model the minimum equations for the phases are also solved by a 
single triangle. When this triangle is imposed on the equations for the moduli it was found 
12] that the solution is always a saddle point: the Hessian matrix has a negative eigenvalue. 



The scalar mass matrix is given by the second derivatives of the potential, and this negative 
eigenvalue is equivalent to the negative mass squared distinctive of false vacua. One may rely 
on radiative corrections to turn the negative mass into positive |y] , but even then the model 
gives a field which is too light to have escaped detection. In consequence, there is no SCPV 
in the singlet model. 



In the 4D model the equations for the phases are solved by three triangles. When 



this solution is imposed on the four minimum equations for the moduli we found []14j that 
they only depend on two combinations of VEVs and, in consequence, there is no solution. A 
solution can be obtained if the mass parameters of the scalar potential are fine tuned and 
two of the four equations become trivial (just like in the model studied in Section 2). But 
even this fine tuned solution is not phenomenologically acceptable, because it is degenerate 
and predicts two massless fields. Radiative corrections would relax the required amount of 
fine tuning and would give mass to all the fields, but still the model has two particles which 
seem to be too light [f|. Therefore, there is no SCPV in the 4D model neither. The Ifi model 
discussed in section II is in some way a particular case of the 4D model, with an analogous 
negative result. 

In this section we explore further extensions of the MSSM. We shall consider the mod- 
els with three pairs of Higgs doublets (6D model), with two pairs of doublets plus one singlet 
(4D1S model) and with two singlets (2S model). When adding singlets, we will not include in 
the superpotential any couplings with dimensions of mass: fi terms for the doublets or linear 
and bilinear terms for the singlets. When such a term does not appear in the superpotential, 
we will not include the corresponding soft SUSY breaking term in the scalar potential nei- 
ther (we assume that the SUSY breaking mechanism respects the discrete symmetries of the 
superpotential). 

The scalar potential for the neutral fields is in each case 
i=i 

+ ( m 2 24 H\H 4 + m 2 6 H\H 6 + m 2 6 HlH 6 + h.c. ) 
+ ( m\ 2 H X H 2 + m 2 Q H\Hq + 777.32 H3H2 + 77734 

+77736 #3#e + m 2 52 H b H 2 + m 2 4 H 5 H 4 + m 2 6 H 5 H 6 + h.c. ) 
+ l(g 2 + g' 2 )[H 2 + Hi + H 2 - H 2 2 - E\ - H 2 6 } 2 , 



;i7) 



V 4D1S = Yl m \ H\Hi + m\&S+ (m 2 l3 HlH 3 + m 2 2A HlH 4 + h.c.) 

i=l 

+( P12 SHiH 2 + flu SH 1 H i + (3 32 SH 3 H 2 + /?34 SH 3 H4 
+^ S 3 + h.c. )+ I a 12 SH 2 | 2 + I a 12 SH X \ 2 + | a 34 SH 4 
+ I a 34 SH 3 I 2 + I a 12 R X R 2 + a 34 # 3 #4 + A SS | 2 
+ \(9 2 + 9 ,2 )[Hl + Hi - H 2 2 - Hi] 2 , 



(1* 



and 



4 1 

V 2S = E m i S l S i + 2 ( & S * S * + ^ S " Sl + ft 5 i + ^ 2 ^ + h-c. 



i=l 



+ | a 3 5 , i5 2 + ^ 5*2 + «i | 2 + | ^ + a 4 S^Si + a 2 S 2 | 2 



(19) 



In the AD1S model we have rotated the doublets in the superpotential so that «i4 = a 32 = 0. 
In the 25" model we shall consider for simplicity only the singlet sector, since this will be 
enough to prove that there is SCPV. 

In the 6D model we search for a complex minimum of type 

(H i ) = ±=v i e* (2 = 1,6), (20) 

with 5\ = 0. The minimum conditions for the phases <5« give five equations. In these equa- 
tions there appear 15 combinations of masses and moduli mfjViVj. Following the procedure 
described in Section 2 we find that the geometrical solution consists of the ten triangles in 
Fig. (2)0. 

A given value of mf-ViVj fixes the multitriangle solution: 

m 2 l2 viv 2 = ai - hi — X\ + yi , ml 3 v iv 3 = a 2 - c 2 - w\ + e 2 , m\ 2 v 3 t> 2 = a 3 , 



m l4 ViVi = x 2 - ci - di + Zi , mf 5 viv 5 = f 2 - b 2 - d 2 + w 2 , m 24 v 2 v 4 = x 3 , 
ml 6 viVQ = fi - y 2 - z 2 + ei , m 2 52 v 5 v 2 = b 3 , ml 4 v 3 v 4 = c 3 , m 2 2G v 2 v & = y 3 , 
m 36^3^6 = e 3 , = w 3 , m 2 m v A v 6 = z 3 , m 2 54 v 5 v 4 = d 3 , ml 6 v 5 VQ = f 3 



(21) 



We can now choose a particular 10-triangle and find the mass parameters which correspond to 
that solution. Once we have these parameters we find the second derivatives of the potential 
to check that the solution is indeed a minimum. It turns out that the minimum is always 
degenerate: any 10-triangle solution occurs for a set of mass parameters giving two massless 
eigenstates (in addition to the goldstone boson of the hypercharge) . The situation here is 
completely analogous to the 4D model. The minimum equations have in general no solution. 
A solution is obtained only for a fine tuned choice of mass parameters that renders two of the 
six moduli equations trivial. Then there appear the distinctive two massless eigenstates. We 



4 In a generic SUSY model with 2n Higgs doublets one can choose (n-l)(2n-l) triangles with all possible 
pairs of phases. These will have n(2n-l) independent sides, equal to the number of mfj that appear in the 
equations. 



have generated numerically random 10-triangle solutions and have obtained always the same 
type of nonacceptable minimum. In consequence, we conclude that there is no SCPV in 4D 
and 6D models: complex minima will require the presence of singlet fields. 

Next we study the 4D1S case. The VEVs can be written 

W = -^V* (i = l,4), 
(S) = v 5 e l5 \ 

(22) 

with 5i = 0. The minimum conditions for the phases <5j define three equations. It is convenient 
to rename (5 2 + £5) — > S 2 , (64 + 65) — > £4, and 35$ — > 5&, and 5 3 remains unchanged. In Fig. 3 
we plot the combination of triangles that solves the 5 equations. The sides in these triangles 
are related to the parameters of the scalar potential: 

Pi2ViV 2 v 5 = a 1 + b 1 +x 2 , p 32 v 3 v 2 v 5 = a 3 , m 2 2A v 2 v 4 = y 3 - b 3 , 
Puv 3 v 4 v 5 = z\ - c 3 , /3i 4 f i^4^5 = c 2 - b 2 , m\ 3 viV3 = o 2 - yi + c x , 
a 12 Xv!V 2 vl = -x 3 , a 34 ;\v 3 V4vl = z 3 , 

^ol\ 2 ol 3 4V X v 2 v 3 V4 = -y 2 , fevf = -xi - z 2 . (23) 

The six triangles in Fig. 3 depend on the ten independent distances above. We proceed like 
in the 6D model, choosing numerically a particular multitriangle solution and adjusting the 
parameters in the potential in order to have that minimum. Then we check if the point 
is really a minimum and we find its properties. We obtain that, although the point given 
by this method has zero first derivatives, it is never a minimum. Like in the singlet model 
analyzed by Romao Jl2|, here the Hessian has the negative eigenvalues that characterize a 
saddle point. For all the (random) cases that we have produced we find negative eigenvalues, 
never a minimum. We conclude that the extension of the MSSM with one singlet plus an 
extra pair of doublets does not offer the possibility of SCPV neither. 

Let us finally consider the 2S model in Eq. (|T9|). We search for complex minima of 

type 

(Si) = Vie* (i = 1,2). (24) 

The two minimum equations for the phases are solved by the set of four triangles in Fig. 4, 
where the six independent distances in the triangles are related to the six independent pa- 
rameters in the scalar potential: 

P3V1V2 = -ci , favivj = -a 2 , 3^it»i = 2c 2 + a 3 + c 2 — - a 3 — , 

Ci a 2 



3p 2 Vn = b 2 - xi + 627^ + xi— , 
03 x 2 

Viv 2 (a 3 a4vl + 2o.\(X' i v\ + a^cn^vl + Z^a^vl) = b 3 , 
2v\v\{(X\Cik + a 2 a 3 ) = —d 2 . 



(25) 



From the triangles we can construct the solution with v\ — 1, v 2 — 1.5, Q\ = n/6 and 
9 2 = 7r/12. This minimum corresponds to f3i = 0.61, (3 2 = 0.05, (3 3 = —0.33, /3 4 = —0.13, 
ai = 2.8, a 2 = 1.3, a 3 = 1, a 4 = —0.58, = —25 and m\ = —6.1. 

The spectrum in the scalar sector is found diagonalizing the 4x4 matrix 



M 



( 


2.8 


-1.0 


-5.8 


4.3 


\ 




-1.0 


0.77 


2.6 


-2.0 






-5.8 


2.6 


91. 


6.7 






4.3 


-2.0 


6.7 


20. 


J 



(26) 



The eigenvalues give m\ = 9.6, m 2 = 4.6, m 3 = 1.1 and m 4 = 0.56. This particular case 
proves that there is SCPV in SUSY models containing two singlet fields. 



4 Conclusions 

It is well known that the standard model and its minimal SUSY extension do not allow for 
SCPV. To know whether in a more complicated Higgs sector CP can be broken spontaneously 
requires solving equations that, in general, are difficult to handle. We have shown that one 
can always build a combination of triangles which is a solution of the minimum equations for 
the complex phases. In simpler cases (3D, 2D IS, ft and 4D models) the triangles are enough 
to solve also analytically the minimum equations for the moduli, and in more complicated 
cases they help to find numerical solutions. 

Using this method we have analyzed the possibility of SCPV in SUSY models. Despite 
the large number of arbitrary parameters present in these models, we find that a Higgs sector 
with only doublets does not provide SCPV: if all the parameters in the Lagrangian are real, 
then the Higgs VEVs can not be complex. This result has been proven for the 4D model, the ft 
model with sneutrino VEVs (analogous to a 5D model), and the 6D model. In consequence, 
SUSY scenarios for SCPV require singlets. We find that, if the singlets do not introduce 
dimensional parameters (i.e., no linear or bilinear terms in the superpotential), one singlet 
is not enough to generate SCPV: the 2D1S and the 4D1S models have always real minima. 



The MSSM extended with two gauge singlets would be the minimal SUSY model where CP 
violation can be generated spontaneously. 

Other possibilities consistent with a soft origin of CP violation would be constrained 
by this result, like the 4D model in Ref. [0]. There all the parameters in the Lagrangian 
are taken real except for the mass terms of the Higgs fields. It is argued that these masses 
could appear at higher energy scales from large VEVs of singlet fields weakly coupled to the 
Higgs doublets, and thus they can in principle be complex. Since our analysis is valid also 
in scenarios with a hierarchy between singlet and doublet VEVs, it follows that the model at 
the large scale must include at least two singlets. Other models that would require at least 
two singlets to obtain soft CP violation within specific supersymmetric models can be found 

in § n, m. 

We would like to emphasize that our method to find complex minima is not restricted 
to SUSY models. Any potential with all the parameters real will give minimum equations for 
the phases of the VEVs that involve only sines (no cosines) of the phases. In this case, using 
the sine law one can define a combination of triangles that solve the equations. This is true 
for all supersymmetric potentials regardless of the CP properties of the fields, as well as for 
all nonsupersymmetric potentials with CP even fields and even for the physically interesting 
minimal extension of the SM with one CP odd real singlet. Of course, there is always the 
trivial solution with all phases equal to zero and no SCPV. But the search for complex minima 
beyond the simplest cases seems almost unworkable unless a method like the one described 
in Section 2 is used. 
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Figure 1: The six triangles that make up the geometrical solution for the minimum phase 
equations for the ft model. The lower right triangle is the solution for the 3 Higgs doublet 
extension of the SM. The relations between the sides of the triangles and the mass parameters 
are given in the text. 




Figure 2: 6D model. 




Figure 3: 4D1S model. 




Figure 4: 2S model. 



